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Abstract 

Holographic superconductors have been studied so far in the absence of dynamical elec- 
rS \ tromagnetic fields, namely in the limit in which they coincide with holographic superfluids. 

C^ ■ It is possible, however, to introduce dynamical gauge fields if a Neumann-type boundary 

condition is imposed on the AdS-boundary. In 3 + 1 dimensions, the dual theory is a 2 + 1 
dimensional CFT whose spectrum contains a massless gauge field, signaling the emergence of 
a gauge symmetry. We study the impact of a dynamical gauge field in vortex configurations 
where it is known to significantly affect the energetics and phase transitions. We calculate 
the critical magnetic fields Hd and Hc2, obtaining that holographic superconductors are of 
Type II {Hci < Hc2)- We extend the study to 4 + 1 dimensions where the gauge field does 
not appear as an emergent phenomenon, but can be introduced, by a proper renormaliza- 
tion, as an external dynamical field. We also compare our predictions with those arising 
from a Ginzburg-Landau theory and identify the generic properties of Abrikosov vortices in 
holographic models. 
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1 Introduction 

The AdS/CFT correspondence has become a powerful tool to study strongly-coupled systems in 
different environments. Very recently its applicability has been also extended to condensed matter 
systems. In Ref. [Tj a gravitational description of a superconductor was proposed whose properties 
have been extensively studied in the last couple of years j2]. In these models a gauge U(l) 
symmetry is broken by a scalar field that turns on near the black-hole horizon. This corresponds, in 
the dual description, to a global U(l) broken by a scalar condensate. Therefore, strictly speaking, 
these models describe either a superfluid [3j or a superconductor in the gauge- less limit. 

The reason for the absence of a dynamical gauge fleld in previously studied holographic su- 
perconductors is the chosen AdS-boundary condition for the U(l) gauge fleld. In most of these 
studies the gauge fleld was chosen to be frozen at the AdS-boundary by imposing a Dirichlet 
boundary condition. One can make, however, the gauge fleld dynamical if one instead imposes a 
Neumann-type boundary condition at the AdS-boundary. In this article we will make use of this 
option to study the role of dynamical gauge flelds in holographic superconductors |1|. 

In a 3 + 1 dimensional AdS space it is known that we can impose either a Dirichlet or a 
Neumann AdS-boundary condition to quantize a gauge theory, being both related by an S-duality 
[1]. In the Neumann case, one flnds a massless gauge fleld in the spectrum that, by means of the 
AdS/CFT correspondence, can be considered to arise from a 2-1-1 dimensional CFT. It is therefore 
an emergent phenomenon. In 4 -(- 1 dimensions, however, a Neumann AdS-boundary condition for 
the gauge fleld is not well-deflned since it leads to a non-flnite Hamiltonian. This will require, as 
we will show, to regularize the theory and absorb the divergencies in local counterterms. In this 
case, the gauge fleld will not be an emergent phenomenon but just an external dynamical gauge 
fleld coupled to a 3 + 1 CFT [5]. 



^ Neumann boundary conditions have been previously considered in holographic superconductors to study 
systems at fixed charge density. In these cases, however, the studied systems are homogeneous and therefore the 
dynamical electric field vanishes. 



An alternative way to understand the distinction between a gauge field in a 2 + 1 and 3 + 1 CFT 
is to look at the zero mode of the Kaluza-Klein expansion of the gauge field in the holographic 
superconductor model at temperatures T bigger than the critical temperature Tc. For a stationary 
vector potential Oj we find, after integrating over the extra dimension, a kinetic term given by 
J d'^xJ-'fj/^iel), where J-^j = diaj — djai and 

forrf = 2 + l, — = -ln(27rT)|^^o for d = 3 + 1 . (1) 



el Ang^T eg g^ 

Here g is the gauge coupling in the AdS model and L is the AdS radius. In 2 + 1 dimensions this 
massless gauge boson mode has a finite norm and therefore remains in the spectrum, while for 
d = 3 + 1 this is a non-normalizable mode and disappears from the set of dynamical degrees of 
freedom. To keep this mode in 3 + 1 dimensions we must then make its norm finite, for example 
by adding local counterterms. 

The impact of a dynamical gauge field in superconductors is expected to be important in inho- 
mogeneous configurations. For this reason we will concentrate here on the vortex configurations 
of the holographic models. We will explicitly analyze the cases d = 2 + 1 and d = 3 + 1, introduc- 
ing, when studying the d = 3 + 1 holographic superconductor, local counterterms to render the 
norm of Oj finite. This will allow us to explicitly see that the dynamical magnetic field B plays 
an important role in reproducing some of the known features of superconductor vortices, such as 
the exponential damping of B inside the superconductor. We will also show that the dynamical 
gauge field changes the vortex configurations of the holographic models, making the energy and, 
correspondingly, the first critical magnetic field Hd independent of the sample size, as expected 
from a true Abrikosov vortex. The properties of the new configurations are qualitatively similar to 
those arising from a Ginzburg-Landau (GL) theory [6], although we find important quantitative 
differences in the size of the vortex core, the profile of B flowing through the vortex and Hd- 

The effect of an external magnetic field on holographic superconductors has been considered 
in [3 [S] and vortex solutions have been studied before for rf = 2 + 1 in Refs. P [TU] . In all these 
cases there was no dynamical electromagnetic (EM) field, and therefore the vortices were not true 
Abrikosov configurations but just superfiuid vortices. Only Ref. [H] showed, for (i = 2 + 1, a non 
trivial profile for the magnetic field of the form of a vortex magnetic tube; it is however unclear 
the origin of this magnetic field and its relation with Abrikosov configurations. 

The organization of the paper is as follows. In Section 2 we describe, similarly in spirit to 
Ref. [12] , effective field theories from which we can obtain model- independent properties of super- 
conductors and superfiuids, and their vortex configurations. This will also help us to make contact 
with the GL predictions q In Section 3 we present the holographic model and explain how to in- 
troduce dynamical gauge fields. We then focus on the holographic superfiuid and superconductor 
(Abrikosov) vortex in rf = 2 + 1 and d = 3 + 1 dimensions, comparing them with those of the 
GL theory. We calculate the energy of these configurations to find the critical magnetic fields Hd 
and Hc2 and show that the holographic superconductors are always of Type II. In Section 4 we 
present a summary of the results and other concluding remarks. 



^What we mean with GL theory in the case of superfiuids is the limit of frozen magnetic fields in the GL theory 
for superconductors. 



2 Effective theories of superfluids and superconductors 

We are interested in the effective theory for time-independent configurations of a U(l) gauge 
field a^ = (ao,ai), where i,j = l,...,d — 1, and a scalar field $ci whose non-zero value will be 
responsible for the U(l) breaking. The effective action at finite temperature T for a^ and the 
order parameter $ci, obtained after integrating out all the other fields of the theory, depends on 
an effective Lagrange density constructed from gauge-invariant operators: 

r = /3 fd''-'xC,s , C,s = Aff(^'j, |A$cir, |$ci|, ...) , (2) 

where Dj = dj — iaj and /3 = 1/T. Eq. (jj]) is defined in some renormalization scheme. We will be 
assuming that this theory depends only on two mass-scales, fi (that later we will associate with a 
chemical potential) and the temperature T. 

The generic effective theory given by Eq. (J2]) simplifies in two limits. In the limit of small fields 
(as compared to fi and T), this theory approximates to the GL theory 



GL 



/3y"ci'^~^^-^^ + |A$GLp + V^GL(|$GL|)}. (3) 



The GL field <I>gl is defined to be canonically normalized, $gl = V^'^ci, where Hq is a positive 
constant, and 

Vgl = -tt^I^glI' + &gl|$gl|'. (4) 

2eGL 

This approximation becomes reliable, for example, close to the critical temperature T < Tc{^) 
where the "condensate" ^ci has a small value. The other useful limit corresponds to slowly varying 
fields, which implies that -Dj$ci and J-'ij are small and 

r ^ /^/^'^'^^d'^ciDJ^^^^j^^^- + |A$ciP + w^(l$ci|)} , (5) 

where h, W and e are generic functions of |$ciP- In the limit of small fields, we obtain the GL 
theory: /i(|$,i|) ^ h{0) = /iq, W{\^,i\) -> VGL(|$GL|)//io and e^d^.i]) ^ 6^(0) = hoel When 
Eq. (IH]) and/or Eq. ([5]) are applicable they can be used to extract model independent features of 
superconductors and superfiuids. 

Consider first the case where we are at large temperatures T > Tc(yu); here the condensate $ci 
is zero, corresponding to the "normal" phase. By decreasing the temperature, T < Tc(/i) at zero 
magnetic field, the modulus of the scalar field ipci = \^c\\ will get a nonzero constant value ipoo- For 
this homogeneous configuration the effective action in Eq. ([5]) is obviously a good approximation 
of the theory Ifl, and the value of -ipoo is determined by the minimum of the potential V = hW: 

dV 

(V^oo) = . (6) 



dip, 



cl 



^Notice however that, generically, this is not the case for the GL theory in Eq. ([3]). 



This configuration corresponds to the superfiuid/superconductor phase. Two important parame- 
ters describing these systems are ^ and A, defined as 

1 1 d^V 1 

These quantities exactly correspond to the inverse mass of the scalar ip^x and Oj respectively. 

In this work we will be considering time-independent vortex configurations with cylindrical 
symmetry as the main example of our theoretical framework. We define (r, 0) as the polar coor- 
dinates restricted to < r < i?, < < 27r. We will always consider the case ^ <^ R and, in the 
superconductor case, also \ <t^ R. We take the Ansatz 

a^ = a^{r) , $,i = e*"Vci(r) , (8) 

where n is an integer, and all other gauge components are set to zero. For n ^ 0, the fields a^{r) 
and ipciir), satisfying the equations of motion from the Lagrangian in ([2]), describe a straight 
vortex line centered at r = 0. If we insert the Ansatz ([H]) into the action in ([S]) we obtain 



r ~ 27rl^'^"3/3 r^drrh{ilj,i)\ ^ J^ . M a^f + {dri^,,f + -An - a^fi^l + W{ijA , (9) 



R 

where V^'^ is the volume of the space orthogonal to the plane {r,(f)). Here the current is given 
by0 

J, = -^^ = 2M^ci)(n - aMi + rdr (^^^^«^) • (10) 

In the vortex case [n ^ 0) ipci goes to ipoo far away from the vortex center; this corresponds to the 
physical fact that a vortex line destroys superfiuidity /superconductivity only in a region close to 
its center. The details of the vortex configurations depend on whether the field a^ is dynamical 
as in the superconductor case, or just a non-dynamical background as in a superfiuid system. We 
consider the two cases in turn. 

2.1 Superfiuid vortex 

For superfiuids the modulus and phase of $ci are respectively associated with the density noand 
velocity Vi of the superfiuid. In the limit of slowly varying fields, Eq. ([5]), we define them as|j 

n,(|<l>,i|) = 2|<l>eip/i(|<l'ci|) , V, = 5.Arg[<l>,i] . (11) 

The field a,^ is not dynamical; it just represents an external angular velocity performed on the 
superfiuid. This is implemented by working in a rotating frame with a constant angular velocity 
fl = a^/r"^. In going from the static to the rotating frame the angular velocity of the superfiuid is 
changed accordingly: v^ ^ v^f, — Qr'^. The current is then given by J^ = ns(f<^ — f2r^). 



'^ Notice that we have defined the current to include the kinetic term of the gauge field. This is done in order to 
facilitate our treatment for both, dynamical and non-dynamical gauge fields. 

^Right dimensions are obtained by putting appropriate powers of the "boson" mass causing the superfluidity. 



Superfluid dynamics coincides with those of a superconductor in the hmit in which the EM field 
is frozen to certain values. This is achieved by taking the limit e — ?• while keeping the external 
magnetic field B = dr-a^/r constant. In this limit the correspondence between the superfluid and 
the superconductor systems is given by 

n^ B/2 , L^^ 2M , (12) 

where L± is the angular momentum and M the magnetization of the system in the direction 
perpendicular to the (r, 0) plane. In the rest of this section we will use the superconductor 
notation. 

Vortices correspond to configurations with n 7^ where ^/'ci varies from zero (at r = 0) to ipoo 
at large r. The exact solution depends on the specific effective action and therefore it is very 
model dependent. Nevertheless, we can obtain the behavior of ipci in the limit r — )■ and for large 
r. Indeed, for r — )■ 0, the condensate goes to zero and the GL action can be applied to obtain 

V'ci oc rl"l . (13) 

For large r, we can use Eq. ([9]) to obtain, in the absence of rotation (5 = 0), 



^cl - ^c 



1 2^^ f., ^00 dh 



(14) 



showing that ~ ^ gives the size of the vortex core radius. For B = the free energy per unit of 
volume V^'^, Fn, is dominated by the third term of Eq. (|9]): 

[^ dr [^ dr 

Jo ''^ Ji "^ 

that depends logarithmically on the size of the superfiuid sample R. This shows that superfiuid 
vortices are not finite-energy configurations in the limit i? — )■ 00. For 5 7^ 0, we have to consider 
the free energy as a function of the angular velocity, obtaining lI 

F^{B) = F„(0) - / Mn{B)dB , (16) 



where Mn{B) is the magnetization (angular momentum from Eq. 0121) ) of the n- vortex configura- 
tion: 

Mn = Ti I drrJ^. (17) 

The value of Mn is approximately given by 

/R / r^ \ f nR'^ R'^ \ 
drr(n-—B\-nnsi^oo)i^ ^5j , (18) 



^In the superfluid case, this is the correct expression for the energy calculated in the co-rotating system with 
respect to the container. 



that leads to 

F^{B) ^ Fo{B) + nn,{^^) (nHn{R/0 - I^R'^b] . (19) 

From this formula we can easily calculate the critical angular velocity B^i above which the vortex 
configuration is energetically favorable. This is given by the B field at which Fi = Fq: 

5,1 :^ J^ In (R/O ■ (20) 

We observe that Bd — )■ when the size of the sample goes to infinity, that is i? — )■ oo. 

By increasing B, more and more vortices are formed up to a critical value Bc2 at which the 
normal phase is favorable. In the limit B — )■ Bc2 the condensate goes to zero, and the GL theory 
can be apphed. One obtains, with a standard textbook derivation, 

1 



Bc2 = 7T^ . (21) 



2.2 Superconductor vortex 

For superconductors, aj, and correspondingly the magnetic field B, are dynamical fields []. Super- 
conductor vortex configurations are therefore described by two fields, ipci and a^. The value of a^ 
varies from zero (at r = 0) to n at infinity, canceling the logarithmic divergence in Eq. ( IT5l) and 
making the vortex energy finite in the limit i? — )■ oo. 

Like in the superfiuid case, we can obtain the behavior of the fields at small and large r in a 
model independent way. At small r, the condensate drops to zero and the GL action can be used; 
in this limit one can derive 

ifjci (X rl"l , a^<xr\ (22) 

At large r the situation is more complicated. If one uses Eq. (jH]) it is possible to show that the 
fields have the following large r behavior 

7/>,i ~ 7/>oo + ^e-"-/^' , a^c^n + aiv^e"''/^' , (23) 



with ^' = ^, A' = A and ipi and ai being constants. Nevertheless, Eq. (125]) shows that higher 
derivatives are not negligible with respect to the first and second derivatives that we included in 
Eq. ([H]). Indeed, we have 

<9^a(^ 1 drtt^ dra^ 

fi'^ ^ (A/i)''-^ jj, ^ fi ' 
where we have assumed, based on dimensional grounds, that the scale /i suppresses the higher- 
dimensional operators, and that A is of order l//i. A similar situation happens for ipci- We are 
therefore led to the conclusion that we cannot neglect higher-derivative terms to describe the large 
r behavior of the fields. Including them, the equations of motion can (formally) be written as 

>l(n)^ci^^(V'ci-^^oo) , Af{D)a^c^^{a^-n) , (25) 



^In this case we call the dynamical magnetic field B, while we keep H for the external magnetic field. 



where Ai and A/" are unknown functions and the box operator acts on ipci and a^ as 

D^ci = -dr {rdri/Jci) , Da^ = rdr f -dra^ j . (26) 

Fortunately, the solutions to the equations above are also of the form of Eq. ( I23l) but with A' and 
^' generically different from A and ^. In other words, the effect of the higher- derivative terms 
is just to change the values of A' and ^'. From this large r behavior we can see that the radius 
size of the vortex core and the radius size of the magnetic tube passing through the vortex (the 
penetration length) are, respectively, characterized by ^' and A'. 

To calculate the external magnetic field H at which the vortex configuration is energetically 
favorable we must obtain the Gibbs free energy. This is given in terms of the free energy F by 

G[Je.t] = F- [d^-'xa^Jl,,, (27) 



where Jl^^ is an external current coupled to the gauge field aj. We can relate Jext to the external 
magnetic field H that it produces, through U 

VxH = elJ,,t. (28) 

Then we end up with the following Gibbs free energy per unit of volume V^"^ of the vortex 
configuration 

Gn[H] = F„ - ^ / rdrdct) BH = F^~ ^H , (29) 



where we have used the magnetic flux condition J rdrdipB = 27m and assumed that H is constant. 
The critical Hd is defined as the value of H at which Gi = Gq that corresponds to 

Hci = ^(Fi - Fo) . (30) 

The exact value of Fi — Fq depends strongly on the model and therefore Hd can only be calculated 
once the model is specified. 

The minimum value of H for which the energetically favorable phase is the normal phase is also, 
as in the superfluid case, Hc2 = 1/(2^ql). The superconductors that have energetically favorable 
vortex solutions, that is ifci < Hc2, are called Type II superconductors, while the others are called 
Type I. When the external field is slightly smaller than Hc2 the condensate has a small value and 
the GL theory can be applied to predict that Type II superconductors present a triangular lattice 
of vortices [13] . Superfiuids can be considered as deep Type II superconductors and therefore they 
also present a triangular lattice of vortices. We will show that holographic superconductors are of 
Type II. 



®In this Maxwell equation of the external field we use eo, defined as the electric charge in the normal phase 
{ipci = 0), to guarantee that when T — >• Tc, the magnetic field B approaches H. 



3 Holographic superfluids and superconductors 

The holographic theory that we want to study is defined [H [14] by a charged scalar \1/ coupled to 
a U(l) gauge field Aa in d + 1 dimensions (a, /3 = 0, 1, ..., d) and an action given by 

S = Jd'+'xV^lj^^iR-A)+^c\ , with £ = -ij-2^--l|Z},v^|2. (31) 

Gj\f is the gravitational Newton constant and the cosmological constant A defines the asymptotic 
AdS radius L via the relation A = —d{d — 1)/L'^; moreover we introduced J-'a/s = daA/s — dj^Aa 
and Da = da — iAa- For simplicity, we have not added any potential for the scalar. We will later 
discuss the implications of including these terms. We will work in the limit Gn -^ and (7 — )■ 
taken such that the gravitational effect of C/ g'^ can be neglected. In this limit the metric is given 
by an AdS-Schwarzschild black hole (BH): 

ds' = ^[-f{z)de + dy']+^^dz\ f{z) = l-{-\ , (32) 

where t is time, z is the holographic direction such that the AdS-boundary occurs at z = 0, while 
the BH horizon is at z = Zh and dy"^ stands for the d — 1 dimensional flat metric. Since we are 
interested in the theory at finite temperature, we will perform the Euclidean continuation with 
compact time it G [0, 1/T] where T = d/{47iZh)- 

3.1 The AdS/CFT correspondence and dynamical gauge fields 

This d + 1 dimensional theory has a dual interpretation in terms of a rf dimensional CFT at 
nonzero temperature. The AdS/CFT dictionary relates the properties of the AdS gravitational 
theory with those of the CFT. In particular, the fields A^ and ^ evaluated on the AdS-boundary 
correspond to fields external to the CFT: 

fl/. = A^\^=o , s = ^1^=0 . (33) 

They are coupled to CFT operators through the interaction terms a^J'^ + sO. The operator J^ 
corresponds to the U(l) current of the CFT theory, while O is a. CFT operator charged under the 
U(l) with Dim[(9] = 3(4) for d = 3(4). Having chosen a nonzero mass for the scalar \1/ in Eq. ([ST]), 
would have corresponded to take another dimensionality for O. We do not expect however any 
important qualitative difference for other choices of the mass. The dual CFT theory, if it exists, is 
supposed to be strongly coupled and the limit (? — )■ in the AdS theory corresponds to be working 
at the planar level in the CFT. 

Integrating over the CFT fields, one can obtain the free energy -F[a^, s] from which the vacuum 
expectation values (VEV) of the CFT operators can be extracted. In the gravity side, F[a^, s] is 
obtained from the d+1 dimensional AdS Euclidean action SEidfi, s] evaluated with all bulk fields 
on-shell restricted to Eq. f[33l) : 

F[a^,s]=TSE[a^„s], (34) 

from which we obtain the VEVs of the currents 

W = ^^'-"^./.|.=o , (O) = ^z'-''D,<iI%=o . (35) 



The matching with the effective theory of Section 2 is straightforward: the gauge field Oj of Eq. ( 133|) 
is identified with that in Eq. (|2]), while ( jj) and (C) of Eq. (!35|) are identified respectively with 
— /3~^5r/(5a* and $ci when renormalized in the same scheme. 

In the AdS/CFT correspondence, the external fields in Eq. (I34p are considered to be frozen 
background fields. This is suited for holographic superfluids, but not for superconductors that 
require the presence of dynamical gauge fields coupled to the CFT. It is easy however to promote 
the external a^ field to a dynamical field. This corresponds to integrating over it in the path 
integral: 



G[s, J, 



ext\ 



Tin / Dae 



-PF[a^,s]+jd'i:L 



y^uu+'^^.Jp 



(36) 



where, for generality, we have added to -F[a^, s] a "bare" kinetic term for a^ (e;, denotes the bare 

electric charge), and have coupled it to a background external current J^^ to define a Gibbs 

energy. Working in the semiclassical approximation |j, Eq. f l36|) leads to the Maxwell equation for 

the gauge field a^\ 

1 



{J^) + ^S,^'^^ + J[ 



ext 



0, 



(37) 



where we have used that ( J'^) = — (5F/(5a^. Let us see how the above procedure can be implemented 
in the gravity side. Using Eq. (1351) . we can write Eq. (137|) as the following AdS-boundary condition: 



'd-3 



9' 



^3-d 



T-, 



^=0 



+ -d.r^ 



z=0 






ext 



0. 



(3^ 



This is a boundary condition of Neumann type that, in order to be consistent with the variational 
principle, requires the AdS model to include the following extra terms on the AdS-boundary: 



(Tx 






4e^ 



piV 



fJ.'^ext 



(39) 



2=0 



Therefore the Gibbs free energy is given by the AdS Euclidean action 5'^; including the additional 
terms Eq. (139|) evaluated on-shell with the bulk fields restricted to the AdS-boundary condition 
Eq. dMD. 

In the particular case oi d = 2 + 1, and in the limit where Cb/g — t- oo (not adding a kinetic 
term for the gauge field on the AdS-boundary), one can show that the theory defined by Eq. (1361) 
preserves conformal symmetry. In this case the original CFT and the a^ can be considered as 
part of a new CFT. Another way to understand this result is given in Ref. P]. There it was 
shown that there are two ways to quantize a gauge field in the four dimensional AdS. We can 
either impose a Dirichlet or a Neumann boundary condition at 2; = 0. Each option is associated 
with a different CFT, S-dual to each other, with different global U(l). While in the first option 
(Dirichlet boundary condition) the gauge field a^ is a background field, in the second one (Neu- 
mann boundary condition) the gauge field is truly dynamical [1] . In this latter case the gauge field 
arises from the CFT as a composite state, as shows the fact that its kinetic term is induced by 
the AdS bulk dynamics. In other words, this local U(l) appears as an emerging phenomenon. As 
emphasized in Ref. [1], this CFT, which includes a dynamical gauge field, has also a global U(l) 



^The semiclassical approximation is valid in the limit g — > and Cb — > 0. 



with an associated conserved current given by J^ = e^'^^'d^aa, and should not be confused with 
the emerging local U(l). Here we also observe that the emergence of the dynamical U(l) can be 
understood without using conformal invariance. Even for a warped space different from AdS, the 
massless zero-mode of a gauge field in 3 + 1 dimensions has finite norm, corresponding then to a 
composite state in the dual d = 2 + 1 theory. This is related to the fact that the gauge interaction 
in d = 2 + 1 is a relevant operator and therefore is dominated by IR physics. It is thus possible to 
send Cb/g to infinity in this case. 

For d = 3 + 1, the situation is different. The current ( j^) contains a logarithmically divergent 
piece given by (in the gauge A^ = 0) 



z 



= -d^J^^f, In z 

2=0 



+ ... . (40) 

z=0 



The appearance of the logarithmic divergence was already expected from the calculation of the 
kinetic term of a, in Eq. ([T]). This can also be understood by looking at the dual CFT interpretation 
of the gravitational theory. Indeed, at short distances (smaller than 1/T) this dual theory is a 
3 + 1 dimensional relativistic theory charged under a U(l). At the quantum level an external a^ 
gauging this U(l) receives corrections to its self-energy that in momentum space go as 

U{p')c:,pHn{pyAl) , (41) 

where p^ is the 4-dimensional momentum of the gauge field and A^ is a momentum cut-off that 
regularizes a logarithmic divergence. Therefore a^ is a state of infinite norm. If our intention is 
to keep the external gauge field in the theory we must renormalize it. A possible way to do so 
is to place a UV-brane at finite z > as in Randall- Sundrum models [15]. Alternatively, we can 
absorb the divergence in the local counterterm of Eq. (139!) . i.e., defining the bare coupling Cb as 

^ = ^ H — :7 In z\z=o + finite terms , (42) 

ei 4 r 

where eo denotes here and thereafter our renormalized (physical) electric charge at the normal 
phase ('i/'ci = 0). Contrary to the d = 2 + 1 case, the presence of the gauge field a^ breaks conformal 
invariance; therefore the gauge field cannot be considered an emerging phenomenon but just a 
new external state coupled to the CFT [5]. The same is true for any d> A. 

We are now ready to study models of holographic superconductors. We are interested in 
vortex configurations where, as we said, the effects of dynamical gauge fields are important. We 
will however present first the holographic superfluid vortex configurations for both d = 2 + 1 
and d = 3 + 1. This will be useful to clarify previous results in the literature [9j, showing that 
these holographic vortices fulfill the expectations of Section 2 for configurations without dynamical 
gauge fields. Then, we will present our main result: the Abrikosov superconducting vortex. 

3.2 The vortex Ansatz 

For both the superfluid and the superconductor we will demand 

s = , ao = fj, , (43) 

10 



where /x is a constant. We fix s = since we are only interested in the case in which the U(l) 
symmetry is broken dynamically by the VEV of O. The constant /i plays the role of a chemical 
potential. As shown in Ref. [H [H], a nonzero /i is necessary in order to induce, at temperatures 
smaller than some critical temperature Tc{ii), a nonzero value for (O) and to have the system in 
a superfiuid/superconductor phase. This critical temperature is given by [U | 



T, ~ 0.03(0.05)/i for d = 3(4) . (44) 

Notice that n ^ breaks the conformal symmetry of the system and, together with T, set the 
scales of the model. 

To obtain vortex solutions we take the Ansatz [9l |TT] 



^ = V'(^,r)e^"^ Ao = Ao{z,r), A^ = A4z,r), (45) 

and the other components of Aa set equal to zero. As in the previous section, n is an integer and 
a vortex corresponds to n ^ 0. We will be working in polar coordinates [dy"^ = dr"^ + r'^dcj)^) for 
(i = 2 + 1 and in cylindrical coordinates {dy^ = dr'^ + r'^dip'^ + dyl) for (i = 3 + 1. The equations 
of motion for the Ansatz ( H5l) are given by 



z 



d-l; 



9A-h-M]+-M-dr^) + [^-'-^^-r^]^ = 



/ A„/A , 1a..a„/a , Mo (^0-^ 



-.d-l J r \ f r 



-'-'^^Ttl^^^^ +-^^ f^^^^ -^^^^7^^' = 0: 



.d-3n f9zAo\ , 1 o / ^ . ^ 2/lo 2 



d.{id)+^dArdrAo)-^r = 0. (46) 



We will impose regularity to our solutions. This requires at z = Zh- 

Zh r r^ 

--d.A, + rdr(-drA,)-^-^^^f^^' = 



Zh \r ) z^ 

^ = 0, (47) 

while at r = we must have 

drA^ = , A^ = 0, 
dri) = forn = , ^j = ioi n ^ . (48) 

For a superfluid, as we explained before, the energy of the vortices is sensitive to the size of the 
sample R. We will therefore limit r < R, where R, as we commented before, is taken much bigger 
than the vortex radius. 
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3.3 Holographic superfluid vortices 

For a vortex superfluid configuration a^ is fixed: 

a^ = A^\,=Q = -Br^ , (49) 

where the constant B represents the external rotation (or, equivalently, the external magnetic field 
for a superconductor in a situation in which the magnetic field can be considered frozen). This 
corresponds to a Dirichlet boundary condition at 2; = 0. 

Also we impose the following boundary conditions at r = i?: 

a^V = , drAQ = , A^ = -BR^ . (50) 

These conditions are consistent with the variational principle which is used to derive the equations 
of motion from the action. The first two conditions represent the physical requirement that, far 
away from the vortex center, the solution should reduce to the superconducting/superfiuid phase, 
which is independent of r, while the third one is a simple option compatible with ( l49l) . 

We have solved numerically Eqs. (H6|) with the boundary conditions Eqs. fHTl) . f l48|) . fj43|) . f l49|) 
and f l50|) by using the COMSOL 3.4 package [16]. In Fig. [T]we present the order parameter and 
the current as functions of r for the n = 1 vortex solution obtained from such numerical analysis. 
Our solutions have the right behavior at r — )■ and r — )■ cxd as predicted in Eqs. f lT5]) and f lT^ 
respectively. We notice however that, unexpectedly, the order parameter (O) develops a small 
bump at around r ~ 12//i, especially for the d = 2 + 1 case. 

It is interesting to know whether our results deviate from those of the simple GL theory. For 
this purpose, we must first specify the input parameters, ^gl and 6gl, of the GL model. We fit 
these two parameters from two predictions of the holographic model: Bc2 and (J^) at large r. The 
value of Bc2 is determined in the holographic model as the value of B at which (O) reduces to 
zero everywhere in space. With this value and Eq. fl2T|) we can obtain ,^gl- From the value of ( J,^) 
at large r as given in Fig. 1 we match with the corresponding current in the GL model, which, for 
B = 0, reads J^^{r — )• 00) = 2n |$gl('^ -^ 00) | ; this allows to obtain 6gl- We find 

^GL ^ 1.1 (0.9)^-1 , boL ^ 3.3/x (12.4) , (51) 

ioT d = 3 (4) at T/Tc = 0.3. Once ^gl and 6gl are determined, we can obtain the prediction of the 
GL model for the condensate and the current as functions of r in the n = 1 vortex configuration. 
We show these in Fig. [H We can appreciate that the holographic vortex differs significantly from 
that of the GL theory. In particular, the radius size of the vortex core in the holographic model 
is considerably bigger than that in the GL theory, namely ^ > ^gl- For T ~ T^, however, the 
holographic model, like any other model of superfluidity, should reduce to the GL theory. We 
have checked numerically that the holographic prediction for (O) and ( J^) approaches that of the 
GL theory when the temperature is very close to T^. 

Next, we calculate the vortex free energy F„ of the d dimensional superfluid. It is then possible 
to verify that B^i behaves as predicted by the effective theory approach, Eq. (!20|) . and also that F„ 
follows, to a very good approximation, Eq. (I19p . Similarly, the results from the effective theory of 

12 



<0> 



0.06 



0.04 



0.02 





10 



15 



20 25 



30 



10 



15 20 25 30 



<J,f,> 



<Jd,> 



0.12 



0.08 



0.04 




0.05 



0.03 



0.01 



20 25 




Figure 1: The modulus of (O) and {j^) (up to a factor L'^^^ j g^) as functions of r from the holographic model in 
the n = 1 superfluid vortex solution for d = 2 + 1 (solid lines on the left) and d = 3 + 1 (solid lines on the right). 
In this plot we chose T jT^ = 0.3 and i? = 0. The dashed lines are the corresponding profiles in the GL model. 
Presented in units of /i = 1 . 

Section 2, can explain the results obtained in Ref. [H]. Indeed, taking the value of ris^^oo) = 0.28y/p 
as in Ref. [9], we obtain, from Eqs. ( IT8l) and ( TT9l) . 



M„ ~ ^Anie^p - O.lR^y^B 



Fi 



Vp 



0.9\n{R/^)-0AR'^B 



(52) 



that agrees PI. as well as Bd in Eq. ( !20|) . with the numerical values obtained in Ref. [9]. 

The value of Bc2 as a function of T, that, as explained before, coincides with that of a su- 
perconductor {Bc2 = Hc2), will be presented in Section 13.41 As discussed in Section 12.21 any 
superfluid can be considered as a deep Type II superconductor and therefore, when B is slightly 
smaller than Bc2, presents a triangular vortex lattice. This property has been checked in Ref. [17] 
for a holographic superfluid for d = 2 + 1. Here we stress that the same remains valid also for 
bigger values of d as it uniquely comes from the fact that, when the condensate is small, the the- 
ory is well approximated by a GL theory. In the next section we will show that our holographic 
superconductor is a Type II superconductor and therefore is also characterized by a triangular 
lattice of vortices for H slightly smaller than Hc2- 



-'^'^Here we point out a missprint in the value of /3„ given in Eq. (24) of Ref. [9 : the correct one is 
0.1(0.2)i?4^. 
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Figure 2: The modulus of (O) (up to a factor L'^~^/g^) and B as functions of r from our holographic model in 
the n = 1 superconductor vortex solution for d = 2 + 1 (solid lines on the left) and d — 3 + 1 (solid lines on the 
right). The dashed lines are the corresponding profiles in the GL theory. Presented in units of /i = 1. 



3.4 Holographic superconductor vortices 

To model an Abrikosov vortex we consider stationary configurations that do not possess a dy- 
namical electric field but only a dynamical magnetic field. Therefore at z = we will impose the 
boundary condition Eq. f H51) for Aq and Eq. (I55|l for Ai that, in polar coordinates, reads 



-d-3 



^3-d 



d.A, 



9' 



+ -T^rdr I -drA, 



2 = 



0, 



where we have taken Jf^ 



ext 



0. At r = i? — 7- oo we impose that 
dr'4^ = , drAo = , A^ = 



2=0 



n . 



(53) 



(54) 



From the set of equations ( l46l) and boundary conditions Eqs. (l47j) . (I48l) . (I43l) . (!53l) and (!54l) . we 
can numerically obtain the superconductor vortex configurations. The profile for the condensate 
iP) and the magnetic field B{r) = drA^\z=o/r are given as functions of r in Fig. [21 We have 
chosen T /T^ = 0.3 and Cb/g — )■ oo for d = 2 + 1, while, for rf = 3 + 1, we have taken e?, to satisfy 
e(^^(T = Tc) ~ 1.7 L/g"^. We observed that the fields have the expected behavior at small and large 
r given by Eqs. fl22l) and fl23l) respectively. Indeed, the vortex profile at large r has changed from 
the behavior of Eq. flT^ to that of Eq. (12^ as expected in an Abrikosov vortex with dynamical 
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Figure 3: A and A' as functions of T from our holographic model for d = 2 + 1 (on the left) and d = 3 + 1 (on the 
right). Presented in units oi jjl— 1. 



EM fields. Similar to the superfiuid case, however, the order parameter (O) shows an unexpected 
slight increase at around r ~ 12//i. In Fig. [3] we show A and A', defined respectively in Eqs. (171) 
and (123|) . as functions of the temperature. For T ^ Tc both quantities diverge as expected, since 
in this limit we have ipoo -^ and therefore A'— J-A— i-oo. AsT— J-O, however, we observe that A 
and A' differ considerably, with A' increasing its value at T/Tc ^ 0.3 — 0.4. A priori, this would 
indicate that the magnetic flux tube becomes broader as T goes to zero, since the penetration 
length A' grows. Nevertheless, we flnd that the situation is more complex; as T — )■ the magnetic 
flux develops two cores, one of size ~ 1/fi while the other ~ A'. This unexpected behavior deserves 
further studies. 

In Fig. [2] we also provide the corresponding curves in the GL theory; the parameters ^gl and 
6gl in the GL potential are flxed as in the superfluid case, Eq. ( 15T1) . while the electric charge cq 
appearing in the GL action is determined by using the second deflnition in Eq. ([7]) applied to the 
GL case, that is Aql = I/Iv^CoI^glI'^ — ^ C)o)|), and by requiring Aql to be equal to the value of 
A' of the holographic superconductor. Again, as in the superfluid case, we observe that the radius 
size of the vortex core is bigger in the holographic model than in the GL theory. As expected, we 
find these differences disappear as T —)• Tc. 

The free energy per unit of volume V^^^ of the vortex configuration is, after taking into account 
the kinetic term in Eq. 036 p . given by 



T 



yd- 



:Sf 



27r / drr 



1 {dra^y 
2el r2 



(55) 



where Se is calculated with the appropriated boundary conditions already stated. Contrary to 
the superfiuid case, we have checked that Fi — Fq is finite for i? — )■ oo thanks to the presence of 
the gauge field. 

To calculate the critical magnetic field ifci we follow Eq. ( l30l) . We find that ifci < Hc2 for any 
real value of cq. This implies that for these holographic superconductors there is always a range of 
values of H for which vortex solutions are energetically favorable; the superconductors are always 
of Type II. In Fig. H] we show Hd and ifc2 as functions of the temperature for the same values 
of Cfc as in Fig. |2l Notice that Hd approaches zero as T — )■ 0. This is due to our normalization 
of H in Eq. (|28|1 that makes Hd oc Cq, which goes to zero as T — )■ 0. We can, however, derive 
Hc2/Hci — >■ oo as T — > independently of such normalization. This is a generic prediction of 
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Figure 4: Hd and Hc2 as functions of T from our holographic model for d = 2 + 1 (solid lines on the left) and 
d = 3 + 1 (solid lines on the right). The dashed lines are the corresponding predictions for Hd from the GL theory. 
Presented in units oi ji — 1. 

superconducting CFT. Finally, we compare our results with those arising from the GL theory of 
superconductors. We observe that Hd deviates from the GL prediction for temperatures smaller 
than Tc- 

From the discussion given in Section 12.21 and the fact that our superconductors are of Type 
II, we know that the energetically favorable configuration when H is slightly smaller than Hc2 is 
a triangular lattice of vortices. 

4 Conclusions 

We have shown how to introduce a dynamical gauge field in holographic superconductors to study 
vortex configurations. In d = 2 + 1 this gauge field is part of the CFT spectrum and therefore 
can be considered to be an emergent phenomenon, instead of an external field. We have shown 
that vortex configurations, in the presence of this gauge field, follow the expected properties of 
finite-energy Abrikosov vortices where the magnetic field drops exponentially at distances larger 
than the vortex core. We have calculated the energy of the vortices and the critical magnetic 
fields Hci and ifc2 that determine the intermediate (Shubnikov) phase. In all cases we have found 
that Hci < Hc2 indicating that holographic superconductors are of Type II. For comparison, we 
have also calculated the vortex configurations in the absence of dynamical fields, corresponding 
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to superfluid vortices, and described their properties. 

The vortex configurations found here differ considerably from those arising from a GL theory. 
In particular, the vortex size comes out to be larger, Hd has a different T dependence, and, more 
importantly, the penetration length of the magnetic field differs significantly as T — )■ 0. 

We have extended the study to d = 3 + 1 where a dynamical gauge field has to be introduced 
by a proper renormalization of the AdS-boundary terms. In this case, the gauge field does not 
respect conformal symmetry and is external to the CFT. In spite of this, the vortex properties are 
found to be very similar to the d = 2 + 1 case. 

Although we have focused on vortex solutions, the method described here to introduce a 
dynamical gauge field is general and can be used in other situations. For example, one could study 
the behavior of the EM field near the surface of a finite size superconductor or in the junction 
between two superconducting samples in the presence of the Josephson effect |18]. Moreover, it 
would be interesting to extend the present analysis to 3+1 dimensional gauge fields sourced by a 
2+1 dimensional CFT; this would allow to study interactions between the EM field and layered 
superconductors. Our studies can also be extended to p-wave holographic superconductors [IS] 
and to holographic models that are dual to non-relativistic scale-invariant theories [20j. 
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